Wrapping Current versus Bulk Integer Quantum Hall Effect in Three Dimensions 
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Surface electron currents are studied for the integer quantum Hall effect in three dimensions(3D) 
proposed previously by Kohmoto et al. and by Koshino et al. We predict the current wraps the 
facets of the sample with its intensity and direction dictated by 3D Chern numbers which are just 
the quantized Hall conductivities in 3D {crxy,Ozx), so a natural connection exists between the surface 
and bulk currents just as in 2D. An experiment to detect the 3D integer quantum Hall effect through 
the wrapping current is proposed. 
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The chiral edge states of 2D electrons in magnetic fields 
have been extensively investigated for the quantum Hall 
effect (Q HE )EJ. Namely, edge states exist for each Landau 
gap in a finite quantum Hall system, and the Hall current 
carried by them is shown to coincide exactly with one 
calculated with the Kubo formulao for the bulk sample, 
which has been interpreted in terms of topological quan.- 
tum numbers characterizing the quantum Hall currentsB. 

While the QHE is usually conceived as specific to two- 
dimensional systems, it is known that integer QHE can 
occur even in three dimensions (3D) if the spectrumiias 
an energy gap and if the Fermi energy lies in the gapoU. 
Although gaps do not usually appear in 3D, two of the 
present authors have shownD that it is possible to have 
a class of energy spectra having a series of gaps (Hofs- 
tadter butterfly) in 3D lattice systems in an appropri- 
ate condition, where we have a 3D-specific IQHE, i.e., 
each of the Hall conductivities <Jxy^ cTzx quantized for each 
gap. That analysis for the 3D QHE is based on the bulk 
description, and a natural question we can now ask is: 
whether and how surface states appear in this 3D QHE. 
In superlattice systems, i.e., stack of 2D systems, the 
surface current has been intensively studiedJj There, one 
discusses a stack of the chiral edge states, called the chi- 
ral sheath current. Surface states are also discussed for 
the field-induced spin density waves (FISDW) in organic 
conductorsO'tj. In those systems, however, the inter-layer 
hoppings are relatively so small that the surface states 
may be understood as weakly coupled 2D edge states. 
By contrast, here we are talking about the 3D-specific 
IQHE with (Jxyi <^zx quantized as 3D Chern (topological) 
numbers. 

So we have studied here the surface states for 3D in 
general and for the 3D butterfly in particular. We show 
that there exists a surface wrapping current that winds 
around the facets of the 3D sample. A hallmark of the 
3D-specific nature appears as the current on each facet 
flowing obliquely to the crystallographic axes, where the 
current direction is given in terms of the 3D Chern num- 
bers. An interesting observation is that the 3D Hall cur- 
rents carried by surface states are exactly the same as 
what is given by the bulk conductivities. 




FIG. 1. (a) A schematic picture of the surface wrapping cur- 
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jn is the surface loop current density. (b)The 



wrapping current calculated for the microscopic model for the 
3D QHE case (eqs.(|ll|) and 

Although we can expect a 3D sample in a magnetic 
field spontaneously carries a surface current in equilib- 
rium as a 2D sample carries an edge current, the situa- 
tion should be more complex, since the surface current 
will have to wind around various faces of the 3D sample. 
However, we can show, from thermodynamics, that the 
surface current may be expressed in a surprisingly simple 
way in terms of the bulk density of states. Let us take a 
cubic, homogeneous sample for simplicity and put it in 
a uniform magnetic field B. We assume that the surface 
current is uniform over each face of the cube and that 
the effect of cube edges can be neglected. Then we can 
characterize the surface currrent fiowing around the faces 
of the cube by regarding it as a bunch of loop-current 
segments as depicted in Fig.|l|(a). Since the current is 
assumed to be uniform on each face, the bunch is char- 
acterized by only two quantities: the vector n, normal 
to each loop, and j, the loop current intensity flowing 
within a unit height measured along n. If we define the 
loop current-density vector as = jn, the magnetiza- 
tion associated with the current is m = yj'°°P, where V 
is the sample volume. Combining with a thermodynamic 
Maxwell's relation, we obtain 
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where p is the density of occupied states per unit volume 
and ^ the chemical potential. 

Now if we move on to the QHE case in the 3D periodic 
system, by assuming that the Fcajwi energy Ep lies in a 
gap of the bulk energy spectrumcla, we have a relation 
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where — e is the charge of an electron and J is a recip- 
rocal vector of the periodic potential. J depends on the 
gap in which Ep lies and gives the quantized 3D Hall 
conductance through 
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where eyfc is the unit antisymmetric tensor. From eqs.(|^) 
and (H), the surface loop current carried by the surface 
states between E and E + is given as 
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From this we will conclude that there is a quantized wrap- 
ping current. Even at this stage, we can see a unique 
behavior: that the direction (in real space) of the surface 
current should always be perpendicular to a fixed recipro- 
cal vector, which in general deviates from the symmetry 
axes of the crystal and also from B, but which is solely 
determined by the gap in which Ep is located. 




FIG. 2. The energy dispersion for the surface states on the 
X faces in the energy gap {M,N) = (3,2). 'L' and 'R' rep- 
resent the surface states on the left (x = 0) and the right 
{x — Lx) faces, respectively. 

Having looked at the thermodynamics, let us move on 
to a microscopic model for the 3D QHE surface currents. 
We take a noninteracting tight-binding electron system 
in a uniform magnetic field B. Schrodinger's equation 
is — ^V'i- Here is the wave func- 

tion at site i, the summation is over nearest- neighbor 
sites, and 6ij — j- A ■ dl is the Peierls phase fac- 
tor arising from the magnetic flux where A is the vec- 
tor potential with V x A = B. We consider a 3D 
simple-cubic lattice with lattice constants a,b,c and B = 
{0, B sm9, B cos6) lying in the yz plane. In the gauge 
A = {0,Bxcos9,~Bxsm9), y,z become cyclic, and we 



have ijjimn = e*'''!''"''+^''^"'^Fi, where Z, m, n are the site in- 
dexes along x,y,z, respectively. Schrodinger's equation 
then reads 
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are the transfer integrals between near- 
est neighbors along x,y,z, respectively, and {Gb,Gc) ~ 
i{B,b,Byc). 

Now we consider that the system is quasi- lEL i.e., 
tx ^ ty,tz, which is required for the 3D butterflyQ. By 
applying the effective-mass approach for the x direction, 
we have 



[E^idx] ~ 2ty cos{GbX + kyb) 

-2tzCos{~GcX + kzc)\F{x) = EF{x), (6) 

where Ex{dx) = {at^dxY ■ The perturbation (oc ty,tz) 
mixes eigenfunctions {F{x) = e^'^"^}, and energy gaps of 
magnitude A = [iyl^^ltzl^ (in a unit of energy tx = 1) 
open at kx = ±-^{MGb + NGc), where M, N are integers. 
In the quasi-lD limit ty,tz — > 0, Schrodinger's equation 
around the gap (M, N) simplifies into 
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where we have decomposed the wave function into the 
left- and right-moving components, F = around 
kx = ifcp = ±i(MGh + NGc), f = Mkyb - Nk^c, and 
V = Tik-p/m* . 

Now we consider the surface states at th&||ga,p (M, N) 
by extending the discussion for the FISDWeIIIJ, which is 
one way of realizing the energy gaps in 3D systems and 
also niatLiematically similar to the 3D QHE considered 
in refE^Eil. We first consider the surfaces normal to the 
conductive axis (x) in a sample with Q < x < Lx- By 
applying the boundary condition F{Q) = to the differ- 
ential eq.(^, we find a surface state that decays into the 
bulk, 
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where k = {A/hv) siii(Mkyb— Nkzc). The corresponding 
eigenenergy measured from the gap center is 



Emn = -~Acos{Mkyb- Nkzc). 



(9) 



The solution with k > corresponds to the left surface 
{x ~ 0) while K < to the right {x ~ Lx). The energy 
dispersion of the surface state must be plotted against 
{ky,kz) as in Fig.||. They oscillate M(7V) times along 
kxiky), where we can see the stripe-like areas alternating 
for the right and left surface states. Since the Brillouin 
zone is topologically a torus, we can define two winding 
numbers for the stripe along toroidal and poloidal direc- 
tions. They are just M and iV, and these two numbers 
are in fact the Chern numbers in 3D, which correspond 
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to UxyjO'zx as shown below. For 2D Hatsugai has shown 
that the edge states whose energy dispersions correspond 
to wiggly hnes against ky between the Landau subbands 
have topological numbers.tl Hence the present result is 
a natural extension to 3D. Although small ty^tz are as- 
sumed above, the 3D Chern numbers, being topological, 
should be constant for larger ty,tz- 

From this microscopic model, we can actually calculate 
the current density on the surface. The expected value 
of the velocity for the electron for g iven {ky-jkz^ is 
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= ±-\/A2-i?2^^(0,M6,-7Vc), (10) 



where + and — corresponds to the left and right sur- 
face states, respectively. We can see that every state in 
the gap (M, N) has a velocity parallel to a single vector 
{Mh, —Nc). The current density on each surface carried 
by the states between E and E + is then expressed as 
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The derivation of the surface currents on the planes 
_L y or i; that contain the conductive axis is slightly dif- 
ferent. We consider a finite sample with < y < Ly,0 < 
z < Lz- To see which states are mixed, we can write the 
perturbational term (the off-diagonal term in eq.(|^)) for 
the gap (M, N) in a second quantized form in k space as 
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where c^(c) are creation (annihilation) opperators. In the 
Wannier representation for y, z, this becomes 

^ — E ^^l:^+MGt+NGc,m-M,n+N^k^,m,n + ii-C., (13) 

kx,m,n 

where \kx,m, n) is a mixed Wannier-Bloch basis localized 
at (y, z) = {mb, nc) and delocalized along x (so we call it 
a 'chain'). Once TC is switched on, the states at {m,n) 
with fcj; > and (m — M, n + N) with fc^, < are mixed, 
and an energy gap appears at kx = ±^{MGb+NGc)- We 
can immediately see that chains lying within M{N) lat- 
tice constants of the faces _L y{z) do not couple to other 
chains, so these states remain gapless, which is exactly 
the origin of the surface states. The expectation value 
of the velocity along x is equal to ±^^{MGb + NGc) 
with opposite directions between the two sides. The cur- 
rent density for each surface for energies between E and 
E + dfj, is 

djy=±-^ f — ,0,Ojd/x iaty = 0,Ly), 



e /N 

dj =T-^ f — ,0,0 (at z = 0,L, 
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From eqs.(pT|),(14), we can see that the currents satisfy 
Kirchhoff's law on each edge of the sample, so that we 
end up with a current sheet that wraps the whole surface 
(Fig.|l|(b)), and the corresponding loop current segment 
defined above is dj'°°P = f (O, f , f-) d/x. We can see 
that this result is consistent with the formula(Q), since 
the reciprocal vector J (eq.(Q)), calculated originally by 
Montambaux and KohmotcB, is {0, -2TrN/b, -2ttM/c) 
for the present case. 




FIG. 3. (a)The experimental setup to detect 3D QHE 
surface currents. The arrows represent the direction of the 
currents, which is opposite to the motion of electrons. (b)The 
"hot line" in the 3D QHE experiment with (c)the correspond- 
ing picture in 2D. 

In the 3D QHE system, how to measure the conduc- 
tivity tensor experimentally is even less trivial than in 
2D. We propose here an experiment to detect 3D QHE 
surface currents in analogy with the 2D Hall bar exper- 
iment. As shown in Fig|| we attach two pairs of elec- 
trodes, (1,2) _L y and (3,4) ± z. The incoming or out- 
going currents, denoted as Ii, ■ ■ ■ , I4 as depicted in the 
figure, are a function of the chemical potentials of the 
electrodes, /ii, • • • , ^4. We can calculate the currents by 
first assuming that all the currents are carried by the sur- 
face states. For the present electrodes, we are left with 
the surface currents on the two x faces (which we call 
right and left surfaces), whose current density is given in 
eq.(0). 

As for the chemical potential, let us consider the situ- 
ation where all /i's are in the gap (M, N) with M, N > 0. 
Then we can see from eq.([Tl|) that the electrons on the 
left surface flow from the electrode 1 or 3 (source) into 2 
or 4 (drain) with a reverse current for the right surface. 
So when we set /ii = ^3 7^ /l(2 = M4, each x face should 
be in equilibrium with the source electrode (^1 for the 
left surface, 112 for the right) if we neglect dissipations 
in the sample (which we will touch upon later). The net 
currents must be conserved, which implies, for rectangu- 
lar surfaces, Ii — l2{= ly) and Iz = Ii{= —Iz), and they 
are readily calculated, via eq.(0), as 
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The Hall conductivity tensor due to the surface conduc- 
tion becomes 
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where we have put fJ-i — fJ.2 — —sVx with Vx being the 
Hall voltage. If we compare this with the expression for 
the bulk conductivity^, we can establish a relationship for 



the Hall conductivities, a. 
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result does not change whether the currents flow in t 
bulk or on the surface. In the usual 2D QHE, Hatsuga: 



shows that cr^™ 



face 



— (Txy by identifing the connection 



between the topological integers for the bulk and the edge 
states. So the discussion here in terms of the 3D topologi- 
cal numbers shows that this property remarkably extends 
to 3D. 

We can give an intuitive way to understand why sur- 
face or bulk does not really matter. Let us start with 2D 
QHE to consider two possible situations (Fig.^) : Case A 
has an electrostatic-potential gradient over the bulk with 
the potential drop, eV, assumed to be the same as the 
chemical potential difference across the two edges. Case 
B has no potential drop in the bulk while the chemical 
potentials at two edges, /ii,/i2, differ. Physically, the sit- 
uation is determined by, e.g., how the local equilibrium 
is achieved by inelastic processes as studied by, Ando in 
a numerical calucation of the potential profilec3. We can 
now question how the a 's in two pictures are connected. 
If we put /^i — /i2 — eV , we can envisage that the two 
cases cross over to each other continuously as in Fig.^. 
The only assumption is that we can neglect the scatter- 
ing across the edge states on the two sides. The current 
in case A can-be divided into the bulk current /^""^ and 
edge currentai3. The edge currents on the left and right 
edges have opposite directions but the same intensity be- 
cause the chemical potential at either edge is assumed to 
have the same energy difference (denoted by /io in the 
figure) from the Landau level center. In case B, the con- 
duction is entirely due to the edge currents /J^'^^" , /j'^^'^ , 
which are different because /ii 7^ /i2- Since the total cur- 
rent must be preserved as we go from case A to B, we 
have /J^'^sc _ j^'^sc _ jbuik^ which is in fact the equality 
in question. If we regard the 3D system as a stack of 
current segments (or 2D Hall "bars") as in Fig.^(b), this 
argument can be extended, since currents across the adja- 
cent "bars" are absent. Then we reproduces the property 

^surface ^bulk 

As a final comment, in the 2D Hall bar geometry it has 
long been recognized that there are two hot spots where 
the chemical potential has to drastically drop from ^1 
to ^2 dissipatively. In our 3D geometry, the hot spots 
should become two "hot lines" as shown in Fig.S(b). 



Given that the required magnetic field (< 40-.I^ is well 
within experimental feasibility as estimated iroEj, exper- 
imental detection of the wrapping currents and hot lines 
should be interesting. 
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Case A 




FIG. 4. An intuitive picture which explains the correspon- 
dence cr'="'^'^'='= and o-'^""^ in the 2D Hall bar. 
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